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Abstract 

The theoretical basis of continuum percolation has changed greatly 
since its beginning as little more than an analogy with lattice systems. 
Nevertheless, there is yet no comprehensive theory of this field. A ba- 
sis for such a theory is provided here with the introduction of the Potts 
fluid, a system of interacting s-state spins which are free to move in 
the continuum. In the s — > 1 limit, the Potts magnetization, suscepti- 
bility and correlation functions are directly related to the percolation 
probability, the mean cluster size and the pair-connectedness, respec- 
tively. Through the Hamiltonian formulation of the Potts fluid, the 
standard methods of statistical mechanics can therefore be used in the 
continuum percolation problem. 
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1 Introduction 



The theoretical basis of continuum percolation [|l|] changed drastically during 
the 1970's. In their seminal 1970 paper, Scher and Zallen || considered 
continuum percolation as an extension of lattice models. Noting that the 
critical volume 4> c is approximately universal, i. e., essentially independent of 
the lattice structure, they suggested it might be carried over into continuum 
systems. 

Such a view of the continuum as the representation of what is lattice- 
independent sits well with the ideas of the renormalization group where the 
continuum often represents merely a scale at which the underlying lattice 
blurs and vanishes. 

Yet continuum percolation turned out to be much richer than suggested 
by this conception. In the first extensive simulations of continuum systems, 
Pike and Seager || found that while the critical volume concept works fairly 
well for spherical objects, it fails completely when the system contains elon- 
gated bodies. This and similar works |§] prompted Balberg et al. || to 
suggest that rather than the critical volume, the proper universal quantity 
might be a critical excluded volume. This concept, however, cannot be re- 
duced to some underlying lattice structure. Instead, it requires considering 
the continuum as fundamental rather than a large scale approximation. 

Fruitful as the excluded volume concept was, it nevertheless had clear 
limitations, since the critical excluded volume is not in fact a truly universal 
quantity. The shape of the objects in the system clearly influenced the critical 
point — i. e., the critical density, p c , at which an infinite cluster first appears 

- in a complicated way. 

In 1977, Haan and Zwanzig || managed to push to its limits the analogy 
between continuum and lattice percolation. Relying heavily on graph theory, 
they managed to expand the mean cluster size, S, in powers of the density. 
From this, they obtained values for the critical density which were in good 
quantitative agreement with the results of simulations. However, this line of 
developments was reaching its limits as it became clear that continuum per- 
colation could be a complex interplay involving not only the shape of objects 

- which determined the binding criterion - but also possible interactions be- 
tween these objects This new element was beyond accommodation into 
the analogy with lattice models. 

The breakthrough which gave continuum percolation an independent the- 
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oretical basis came in two papers by Coniglio et al. M. Building on previous 
work by Hill || , these researchers showed that the mean cluster size could be 
related to a pair-connectedness function which had an expansion in powers of 
the density. This expansion included naturally the interactions. Originally, 
Coniglio et al. developed this expansion in the context of physical clustering 



in a gas, but it was soon extended by analogy to other systems pQ, 11]. Al- 



though the results were mainly qualitative at first, some recent calculations 
have obtained quantitatively good agreement with simulations |L2] . 



In spite of these successes, the theoretical basis of continuum percolation 
is still incomplete. The theory of Coniglio et al. suffers from two essential 
problems. The first is formal. The fundamental quantity in this theory 
is the pair-connectedness function, yet the only mathematical definition of 
it is in terms of a power expansion. This is an unsatisfactory situation. 
Some non-perturbative scheme must underly such a result, but it is absent 
from the work of Coniglio et al.. Another unsatisfactory aspect is that the 
extension of this work to various percolation systems has to be done by 
analogy and reasonable assumption. While this is not a serious problem, one 
would like to see a more formal basis for such uses. The second main problem 
is that the theory of Coniglio et al. does not say anything about the order 
parameter, the percolation probability. This problem is related to the first 
because perturbative expansions are limited to the region of densities below 
the critical point, while the order parameter is non-trivial only at densities 
above this point. Furthermore, the percolation probability is not obtainable 
from the pair-connectedness, and requires therefore a new approach. 

Such limitations also prevent one from approaching some fundamental 
questions. First among these, probably, is related to the continuum perco- 
lation universality class. Computer simulations suggest that continuum per- 
colation belongs to the universality class of lattice percolation. One would 
like to understand this from the theoretical point of view, as well as to as- 
sess whether this is true for all possible interactions and binding criteria. 
Computer simulations are inadequate to address this problem since they are 
greatly complicated by the presence of interactions. Furthermore, obtaining 
critical exponents from a simulation is a notoriously delicate procedure, and 
one cannot claim yet that the universality class of continuum percolation is 
known with certainty. 

The present series of papers is meant to be the basis for a more compre- 
hensive theory of continuum percolation, which is based on an extension of 
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the classic mapping between the Potts model and lattice percolation which 



was invented by Kasteleyn and Fortuin [13]. The present extension includes 



an off-lattice version of the Potts model which I have called the Potts fluid. 
The s-state Potts fluid is introduced in section 2. I show then that, as in the 
lattice case, all statistical averages in the percolation model may be expressed 
as averages in the Potts fluid. In section 3, this mapping is applied to the 
Potts magnetization and susceptibility. The limit s — *■ 1 of these quantities 
yields the percolation probability, P(p), and the mean cluster size, S, respec- 
tively. Section 4 is concerned with the Potts correlation functions. The pair- 
connectedness is shown to be directly related to the Potts pair-correlation 
function. Section 5 sums up the results. 



2 General formalism 

By "continuum percolation" I always mean a system of particles interacting 
through a pair potential v (f, fj) and obeying classical statistical mechanics. 
In addition to their interaction, the particles possess another property, "the 
connectivity", determined by a probability function p(fi,fj) which is the 
probability that two particles located at f and fj are bound, or connected, 
to each other. A cluster is a group of such connected particles. Let us also 
define q(fi,fj) = 1 — p(fi,fj) as the complementary function. 

Note that we do not assume any "locking" mechanism: bound particles 
do not remain glued together thereafter. Indeed, the system does not evolve 
with time at all. The point of view adopted here is that of equilibrium 
statistical mechanics, where all properties are derived from an ensemble of 
"snapshot" configurations. As a result, the properties of connectivity and 
interaction are arranged hierarchically. The interaction alone determines the 
configuration of the particles, which in turn helps determine the connectivity, 
but the connectivity does not in its turn influence the configuration. 

Often, as in systems of permeable objects ||, p(f, fj) takes only two val- 
ues, or 1, but no such restrictions are assumed here. Thus, the connectivity 
state need not be uniquely determined by the geometrical configuration of 
the particles. Only its probability distribution is fully determined through 
the function p(fi, fj). 

The formalism presented here is an extension of the Kasteleyn-Fortuin 
mapping between lattice percolation and the Potts model [|13[]. For the con- 
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tinuum case I will define an extension of the Potts model hereafter called the 
Potts fluid. 

The s-state Potts fluid is a system of N "spins" {Aj}^ each having s 
possible states, and obeying classical statistical mechanics. Each spin has a 
position fi in the continuum and the spins interact with each other through 
a spin-dependent pair potential V(fi, Xf, fj, Xj), such that 



V(fi, X^fj, Xj) = V(i,j) 



U(r i ,f j ) if Xi 
W{f h fj) if Xi ^ Xj 



A, 



(2.1) 



where U and W are arbitrary functions. 

The spins may also couple to an external field h(r) which tends to align 
them in some given state. If we arbitrarily denote this state as "1" , then 
the interaction Hamiltonian is 



N 



Him = ~J2 ^{'Xi)Wi) 



(2.2) 



»=i 



(2-3) 



where 

^ (A) = { -1 if A/l 

For conciseness, two spins in the same state will be said to be parallel 
to each other (though the state need not actually correspond to any spatial 
directions). Two spins in different states will be called non-parallel. 

As in all classical systems [fL4|, the dependence on the momenta can 
be factored out so that all statistical averages depend on the configuration 
integral 



1 

ivT 



E 

{A m } 



dr\ • • ■ dr^ exp 



N 



i>j i=l 



(2.4) 
1/kT is the 



where the sum J2{\ m } is over all spin configurations, and (5 
inverse temperature as usual. 

Now, any continuum percolation model defined by v(i,j) and p(i,j) can 
be mapped onto an appropriate Potts fluid model with a pair-spin interaction 
defined by 



U{i,j) 
exp [-/3W(i,j)) 



q(i,j)exp [-(3v(i,j)} 



(2.5) 



5 



This mapping relates statistical averages of the Potts fluid to statistical 
averages of the percolation model. The most fundamental relation is ob- 
tained for the Potts fluid configuration integral Eq. (|2.4| ). However, because 
the following derivation can be a little confusing, let us anticipate the final 



result. The formal mapping Eq. (2.5) induces a geometrical mapping be- 



tween spin configurations and connectivity states. To every particle in the 
percolation system we assign a spin in such a way that if two objects belong 
to the same cluster in the percolation model, they are assigned the same spin 
(i.e, their spins are parallel) in the corresponding Potts fluid spin configura- 
tion. However, the actual value of this common spin is selected at random. 
Therefore the mapping from the percolation model to the Potts fluid is one 
to many. Any connectivity state corresponds to several spin configurations, 
which differ from one another by the actual value of the spin assigned in 
common to all the particles in a given cluster. Because of the one to many 
character of this mapping, the status of two parallel spins differs radically 
from the status of two non-parallel spins. If two spins are non-parallel, they 
must belong to different clusters in the percolation picture. However, if two 
spins are parallel no conclusion may be drawn as to their connectivity state. 
They might be parallel because they belong to the same cluster. But they 
might also belong to different clusters and have been assigned the same spin 
by chance. Since the common spin of every particle in a cluster is assigned 
randomly, it is quite probable that two unrelated clusters end up with the 
same spin. This difference makes "disconnectedness" a more basic property 
than connectedness in the sense that given a spin configuration, one can 
never deduce with certainty that two particles belong to the same cluster in 
the corresponding percolation model. However, one may be able to deduce 
that two particle do not belong to the same cluster. 

Let us now see how this geometrical mapping arises naturally out of the 
formal mapping Eq. ( |2.5|) . For simplicity, let us assume temporarily that the 
external field vanishes, i. e., h(r) = 0. Introduce the notation 



1 

ivT 



dn ■ ■■df N Q(fi, ...,r N ) 



<3(Fi, ...,fjy) = ex P 

{A m } 



i>j 



(2.6) 

In the expression for Q(l, . . . , N), let us separate all possible configurations 
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A m into those where Ai = A2 and the rest 



Q = e -/»W) J2 ex P 

/ A m \ 

-/w(i,2) J- exp 

/ A m \ 



+ e 



-0 E nu) 

(<J)-(1,2) 

E 

(i>i)-(l,2) 



(2.7) 



where S(ij)-(i,2) means a summation over all pairs except the pair 

(1, 2). In the second term on the l.h.s of Eq. ( |2.7D we can now rewrite 



E ex P 

/ A m \ 



-0 E 

(ij)-(1.2) 



E ex p 

{A m } 



E v(i,j) 

(<J)-(1,2) 



E ex p 

/ A m \ 
\ Al=A 2 / 



-0 E 

(ij)-(l,2) 



(2.8) 



where the sum S{A m } is now performed over all spin configurations without 
constraints. Eq. (|2.7p now becomes 



Q 



e -/3C/(l,2) _ e PW{l,2) 



E ex p 

/ A m \ 
\Ai=A 2 J 



-/3 E ^>i) 

(ij)-(l,2) 



+ e 



-/3W(1,2) 



E eX P 



{A m } 



E V(i,j) 

(W)-(1,2) 



(2.9) 



The mapping Eq. (|2.5| ) then implies that 

Q = p(l,2)e- /3l ' (1 ' 2) E ex P 

/ A m \ 
\ A i=A 2 J 



(i,i)~(l,2) 



+ g(l,2)e- /3 ^ 1 ' 2 ) E ex P 

{A m } 



(i,i)-(l,2) 



(2.10) 
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where we used the fact that p(l,2) = 1 — g(l,2). Repeating this procedure 
for the pair (1, 3), we obtain 



Q = p(l,2)p(l,3)e 



-f3[v(l,2)+v(l,3)} 



E ex p 

\ A 1= A 2 = A 3 J 



-P E v(ij) 

(i,j)-(l,2)-(l,3) 



+ p(l,2)g(l,3)e- /3 [ 1 '( 1 ' 2 ^ 1 ' 3 )] E ex P 

/ A m \ 
\ A 1 = A 2 J 

+ g(l,2)p(l,3)e-^ 1 ' 2 ^ 1 ' 3 )] E exp 

/ A m \ 
\ A 1= A 3 / 



-0 E *^') 

(i,j)-(l,2)-(l,3) 

E 

(i,j)-(l,2)-(l,3) 



+ g(l,2)g(l,3)e-^ 1 ' 2 ) + ^ E ex P 

{A m } 



-/? E 

(iJ)-(l,2) 



(2.11) 



One sees therefore how the geometric mapping comes about. Every sum 
is performed over a subset of spin configurations determined by a set of 
constraints of the form Ai = A2, . . .. All constraints are equalities between 
spin pairs. There are no constraints forcing two spins to be non-parallel. 
Every constraint Aj = Xj also brings forth a factor p(i,j), the probability 
for the pair to be bound. Since a bound pair belongs by definition to the 
same cluster, we see the beginning of the geometrical mapping in which spins 
in a single cluster have parallel spins. To see the full mapping, we need to 
repeat the procedure outlined in Eqs. ( |2T7j)-( |2.lT| ) for all pairs of spins. Let 
us consider one of the sums into which the function Q has been decomposed 
a step before, and consider a pair (m, n). Two possibilities arise: 

1. Previous constraints already determine that A m = A n (for example, 
there could be some k for which A m = A& and A n = A&). Then, neces- 
sarily, V(m,n) = U(m,n) = v(m,n). Hence 



E 



exp 



-P E v&j) 



previous 
constraints 
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-/3v{m,n) 



E 



cxp 



-P E v(ij) 

(rn,n) 



previous 
constraints 



p(m, n) e 



-j3v(m,n) 



E ex p 



prev. 
constr. 



+ q(m, n) e~^ M £ exp 



prev. 
constr. 



-P E v(ij) 

(m,n) 



-P E v(ij) 

(m,n) 



(2.12) 



where we have used p(m, n) + q(m, n) = 1. 



2. Previous constraints do not determine that \ m = \ n . Then the sit- 
uation is as it was for the pair (1,2), and the sum will split in the 
following way 



E ex P 

/ Am \ 

\ constr. j 



-P E v(ij) 

. (ij)— ■ 



= p(m, n) e /M"»>«) ^ e xp 

J Am:- ) 

+ q(m,n)e-^ m ' n ^ ]T exp 

{A m :...} 



-P E 

(rn,n) 



-P E ^(U) 

( m >«) 



(2.13) 



From the point of view of the geometrical mapping, case (1) means that 
m and n already belong to the same cluster by virtue of some other particles 
which link them (perhaps indirectly). Hence, it doesn't matter whether they 
are also directly bound to each other [a case which contributes the factor 
p(m, n)\ or not [ which contributes the factor q(m,n)}. As a result, 

both contributions appear without the addition of any constraint. In case 



9 



(2) on the other hand, the appearance of the factor p(m, n) also requires a 
new constraint which now implies that m and n are bound and therefore 
belong to the same cluster. The factor q(m,n) requires no such constraint, 
since if the particles belong to different clusters, their spin attribution is 
random. 

When all pairs have been covered, the set of constraints of a particular 
sum specifies exactly which particles belongs to which clusters in the original 
percolation model configuration. Since the expression for Q contains sums 
over all possible constraints, it can be rewritten as a sum over all possible 
clusterings of the original percolation configuration. Thus, let us define 

i>j 



(2.14) 

Q= E E P(conn.) (2.15) 

all possible 1 f * m consistent ' 
connectivity \ J witn the 
states J I connectivity 
[ state 

where the sum over all spins is consistent with the connectivity state in the 
sense of the geometrical mapping, i. e., that all particles within a single 
cluster must be assigned the same spin (note, however, that not all particles 
in a cluster need be directly bound to each other). 

From the probabilistic interpretation of the functions p(i,j), q{i,j) and of 
the usual canonical Gibbs distribution, we now see that up to a normalization, 
Eq. (|2.14j) means that 

P( conn.) oc Probability density of finding a configuration^, . . . , r^) 
such that all pairs are bound and all pairs (m,n) 
are unbound. (2-16) 

Returning now to Eq. ( |2.4j ), and restoring the external field h(r), we finally 
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P( conn. 'I 



n p(hj) 

all 1 

bound 
pairs 



n 



q{m,n) exp 



all 

unbound 
pairs 

(m,n) 



Then we can write Q as 



find 



l^Jdl-.-dN £ P(conn.) £ 



exp 



{states} ^ Am ' °^ 



' N 
. i=l 



(2.17) 

Here the symbol Z){>- cl} means a sum over an s Pi n configurations consistent 
with the clustering, i. e., with the connectivity state in the sense of the 
geometrical mapping. 

This expression has the form of an average over all percolation config- 
uration (up to nonimportant normalization factors). It is the fundamental 
relation between statistical averages calculated in the percolation model and 
quantities calculated in the corresponding Potts fluid. As we shall see in a 
moment, it allows us to obtain percolation- related quantities by calculating 
properties of the Potts fluid defined in Eq. (|2.5|) . However, since there is no 
equivalent in the percolation model to the field h, we need to set this field to 
zero. In this case, Eq. ( [2.171) becomes 

Z = — f dl---dN Yl ^(conn.) s Nc (2.18) 

N - / conn. \ 

\ states / 

where N c is the total number of clusters in the relevant connectivity state. 
The term s Nc is then the number of possible assignments of spins to these 
N c clusters in accordance with the geometrical mapping. 

Finally we shall have to get rid of the free parameter s. As in the usual 
Kasteleyn-Fortuin mapping |13|], the interesting case is the limit s — > 1. The 
resulting expression 

Zp = ^- f dl---dN Y, ^(conn.) (2.19) 

iV • J f conn. \ 

\ states; 

is just the normalization factor required for averages performed within the 
percolation model. Hence, 



z = z p l Yl ex P 

\{A m : cl} 



N 



i=l 



(2.20) 
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where the symbol ( ) means a canonical average performed in the percolation 
model, i. e., 

(if (1, . . . , N)) p = -J- / dl ■ ■ ■ dN £ P(conn.) F(l, . . . , N) (2.21) 
iy -^p J /conn. \ 

\ states ; 

where F is some quantity defined in the percolation model. 

Similarly, we need averages performed in the Potts fluid, which will be 
denoted by ( ) s . Given a quantity G defined in the Potts fluid, we have 

(G(1,X 1 ;...;N,X N )) S = ^ J dl ■ ■ ■ dN £ G(l, X 1 ; . . . ; N, X N ) 



x exp 



{A m } 

N 

i>j i=l 



(2.22) 



Eq. (|2.15| ) now implies that in general, for any quantity G(l, Ai; . . . ; N, Ajv) 



defined in the Potts fluid system, we have that 

(G(l,Xi,...;N,X N )) s = / £ G(l,X 1 ;...;N,X N )eK h ^A 

HA™: cl} I p 

(2.23) 

This fundamental relation allows us to translate every quantity in the Potts 
fluid into some quantity defined in the continuum percolation model. 

We now turn to using Eq. ( |2.20| ) to show that the average magnetization 
in the Potts fluid is directly related to the order parameter of the percolation 
model, the percolation probability. 

3 The Potts magnetization and the percola- 
tion probability 

The magnetization of the Potts fluid is defined as in the usual Potts model 
by 

1 d\nZ , . 
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where h is the now constant external field. From Eqs. ( |2.20| ) and ( |2.23| ), we 
have 



N 



M 



i ( N r N 

^TV E E^)exp P^hftMXi 
' {A m : cl} Li=l L , i 



(3.2) 



In order to calculate this expression, we need to characterize more for- 
mally the sum over all spin states. In a given connectivity state, let us 
denote by N r (I < r < N) the number of clusters containing exactly r par- 
ticles (possibly N r = 0). Also, because of their particular role, let us denote 
separately the number of spanning clusters by N s (as we take iV — > oo, N s is 
either or 1). Any connectivity state thus corresponds to a set of numbers 
(Ni, N 2 , . . . , N N , N s ) (other connectivity states may also correspond to the 
same set). For a given spin configuration consistent with this clustering, let 
us denote by k r the number of r-clusters which have been assigned the spin 
1 (k s will represent the spanning clusters assigned the spin 1). Every spin 
configuration generates therefore a set (ki, k 2 , ■ ■ ■ , k N , k s ). Then, the total 
number of spin configurations consistent with the given clustering is 



Ni N 2 N N N s 

EE---E E 

fc 1= 0fe 2 =0 k N =0k s =0 





1) 



iVi-fci 



(s-1) 



N B -k s 



(3.3) 



The factor (s — l^ Nr ~ kr is the number of possible spin assignments to the 
(iV r — k r ) r-clusters which have a spin different from 1. 

Consider now the expression J2iLiip(^i)- From the definition of ^(A), 
Eq. ( [2.3|) , every j-cluster with a spin 1 contributes a term j(s — 1) to this 
expression. Every j-cluster with a spin different from 1 contributes j(— 1) to 
it. Therefore, for a given spin configuration, 



N 

Ev>(a, 

1=1 



(s-1) 



N 



E/ foj n s^i 



N 



X./i.V ; - kj) ■ //.( .v. • • /,-. 

3=1 



N s 

E^'% 



N 



(3.4) 
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where n s is the number of spins in the spanning cluster (we anticipate the 
thermodynamic limit and assume that the spanning cluster, if it exists, is 
unique), and where we have used the identity J2f=ijkj = N. The notation 



J2f=ijkj is shorthand for n s k s + Z/j=i jkj. 
We now have that 



N(s-l) S I 
v ' {\ m : cl} k 



N 



i=l 



N 



E^) 



i=i 



-f3hN 



N a ( Nt 



X 



N(s - 1) £j 



E (^) - * 



(3.5) 



JV S N N 

where II = II • • • x (case I = N s ). 
i=i i=i 

It is easy to prove (e. g., by induction) the two identities 



ki=o V JV/ / 



e Ws + ( s _ i) 



TV, 



fc,=o V v '> 



e^ fts + (s - 1) 



A'; 



A, 



e^ s + (s - 1) 



(3.6) 

iV,e i/%s (3.7) 



Therefore, 

E 

.fc(=0 

and 



n 

/=i 



sift- 



AT-s 

AT) = —N Yl [e ll3hs + (s - 1) 
z=i 



(3.8) 



iV s JV, /. - 



ki 
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x E 



A s 



- n [e^ s + (a - 1) 



A, 



(3.9) 



Combining these two equations, we finally obtain 



E | ex P 



{A m : cl} 



A/ 



A 



E<W 



j i=i 



^ r lA/ AT, 

-,7,!V -Q |" e i/9ft. + (s _ ^ i -AT + E SjW. 



a j(3hs 



1=1 



{ J 3 e^ hs + (s- 1) 



(3.10) 



Recalling that iV = J2j=ijNj and rearranging terms, we have that 



N(s- 



7 {A m : cl} L 



exp 



A 



ph524>(\ 



i=l 



A? 



J i=l 

A s 



.v n[«*+(.-i)r|Dv +( .-i) 



e 

at, ^ j iV} 



2=1 

1 f fs - l) e - s/m 



.7 = 



2=1 



n[ 

Hence 

M=(lt[l + (a-l) e 



E 

3=1 



iV 



1 + (s - l)e-i^ s 



(3.11) 



■ AT S 



-s/32ft 



U=l 



1-e 



-j/3/is 



1 + (s - l)e"J 



(3.12) 



To obtain a quantity directly related to the percolation model, we need 
to take the limit s — *■ 1, which yields 



A. 



A., 



^D = (^E^-^E 



3=1 



iV 



3=1 



A., 



Ei^- 



-3"y9A 



\ 



(3.13) 



Finally, we need to take the thermodynamic limit and set the external 
field to 0. As always, we must be careful in the order of these limits A 
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broken symmetry state will only be obtained if the field remains finite while 
N — > oo. In the limit N — > oo, the spanning cluster behaves differently from 
the finite clusters. Eq. ( |3.13| ) can be rewritten as 

M = 1 - 1 (E jiV.e-^ + V^le-^X (3 u) 

In the limit AT — > oo, the factors exp(—j/3h) remain finite, but exp(—/3hn s ) — > 
0, because the spanning cluster then becomes infinite. It is this property 
which distinguishes it from all finite clusters. Hence 

M=1 -Ji(^^) (3 - 15) 

Finally, we can set h = 0. At this stage, the limits commute and we obtain 
M=l - lim /-YjN,S =1- lim = lim /^\ (3.16) 

J p 

This is exactly the probability that a particle picked at random belongs to 
the infinite cluster, which is by definition the percolation probability P{p)- 
Therefore, 

lim lim lim M = Pip) (3.17) 

h— >0 iV-+oo s— >1 

and the percolation order parameter is directly calculable from the Potts 
magnetization. 

We can calculate the Potts susceptibility most easily from Eq. (|3.15|) , 
before we set the field to (it is easy to see that the derivative d/dh commutes 
with all the operations performed up to that point). Thus, 

x = ^ = (3 lim /^Ej 2 A^A (3.18) 
dh N ^\Nj^ 

3 p 

Taking now the limit h — > and using again that N = jNj gives 

i=i 

TV 

N E jNj 1 
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For densities lower than the critical density, there is no spanning cluster (in 
the thermodynamic limit), so that 




X = /?Jim ) =PS (p< Pc ) (3.20) 



where S is the average mean cluster size, by definition fL6| . 

Therefore, we can calculate the percolation probability by calculating 
instead the magnetization of a corresponding Potts fluid, and the average 
mean cluster size by calculating the Potts susceptibility. 



4 Correlation functions 



The Potts magnetization and susceptibility are directly related to the n- 
density functions of the Potts fluid, which are defined in analogy to the 
n-density functions of a classical liquid [M as 



p (n) (fi,Ai;f 2 ,A 2 ;...,f n ,A n ) = 

-r / dr n+1 ■ ■ ■ dr N exp 

— n)l J 



Z(N 



The normalization is chosen so that 



N 



-0£V(i,i)-/35>(*)tK*i) 

i>j 1=1 



(4.1) 



£ Jdl---dnp^ = s N 



{Ai,...,A„} 



(N -n)\ 



(4.2) 



The n-density functions can also be expressed as canonical averages, 



N 



p {2) (x,a;y } 7]) 



i=l / s 

N N \ 



(4.3) 
(4.4) 



\i=l J=l 
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Now, we can rewrite 



^(A.) = (a - l)5 XiA - £ 5 XuX 



(4.5) 



Using the identity / dx 5{fi — x) — 1, we now have 

= (i s ~ !) £ J dx8{r i -x)8 Xi ,i 
~ \Y1 J dx5(fi - x)J2$> 



I n 
\i=l 



dx 



(s-l)p«(£,a=l)-£p«(a;,a) 



(4.6) 



Since the field h preserves the symmetry between all the spins A ^ 1, p^'(x, a) 
must be the same for all a ^ 1. Hence, 



, N \ | 



M 



iV(s 



i=i 



- / dx 



p {1 \x,l) -p [1) (x,a) (4.7) 



where a denotes some (arbitrary) spin value other than 1. From the normal- 
ization, Eq. ( |4.2| ),it is now obvious that 



J dx p^\x, 1) = ni = number of spins in state 1 

J dx p^(x, a) = n = number of spins in any state 1 (4.8) 



Hence 



M = — (n x - n) 



(4.9) 



is just the excess density of spins in the state 1 over the density of spins in 
any other state. 

From Eq. Q3.2|) , and the definition x — dM/dh, the susceptibility is 



X 



P 



N(s 



N(s 



8 N 



(4.10) 



1=1 
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Repeating the steps leading to Eq. (|4.6|) , we have 



3 N 

E^ 







N(s-l) 



P {1) (x,a) 



(4.11) 



where the sum over a runs over all the s possible spin states. 
On the other hand, 

E iKAiMAi) = E E V>(«M7) E <Wa„ 7 



(4.12) 



where, again, the sum over a and 7 runs over the s spin states. Hence, by a 
derivation similar to the one used for Eq. ( |4.6| ), we have 

(5>(W(Ai)\ 



E E ^("Mt) Y.J d % d y - x) 6(fj - y) S Xi , a 5 Xjn 



a 7 



dx dy 



E E ^OM7) p (2) (z, «; f, 7) 



o 7 



(4.13) 



Therefore 







+ 



N(s - 1) 





(if dy 



E E ^( a M7) p (2) (^ £ 7) 



a 7 



N(s 



— — J dx^2ip 2 (a) p [1) {x,a) 



(4.14) 



let us now relate this expression to the percolation picture, i. e., to the 
limit s — *■ 1, when x — > /35 1 for densities p < p c . For this range of den- 
sities, the symmetry of the system is completely unbroken and p^ 1 ' (x, a) is 
independent of a. From the normalization Eq. (4.2), it follows that 

( dxp {1 \x,a) = - ( df^p (1) (f,a) = N (4.15) 
J s J a 

Now E a ^ 2 {a) = (s - l) 2 + {s- 1)(-1) 2 , so that 



(4.16) 
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and therefore 

3 TT / dxJ2^ 2 (a)p (1) (x,a) = (3 (4.17) 



N(s 

Hence, when s — *■ 1, 



^ ^ £ = 1 + lim — ; / dx dy 

(3 N(s -1)J y 



E E ^( a )^(i)p (2) a; 7) 



a 7 



(4.18) 



Let us now recalculate this quantity by using the general connection be- 
tween averages in the Potts fluid and in the percolation model, Eq. ( |2.23|) . 
In particular, we have 

(^g^(A^(A,)) r ^g( { E^(A.)^)> p (4,9) 

where we have already set h = 0, which is allowed since in the range p < p c , 
the symmetry is unbroken anyway. Let us now calculate the average on the 
r.h.s separately when Aj and Xj belong to the same cluster and when they do 
not. To this end, define a function ft(i,j) as 

, N f 1 if 2,7 belong to the same cluster ,.„„> 
"(•''H otherwise ^ 



Then, from Eq. ( |2.21| ) 



E ^(W(Ai)) 



= -L-[di"-dN E P(conn.) £ ^(A#(A^( M 1 

^•^P- 7 /Conn. \ {A m :c/.} 

\ states j 

+ -±-fdl...dN E ^(conn.) E ^(W(A,) [1 - fi(z, j)] 

ly - Zj P J /COnn. \ {A m :c/.} 

\ states j 

(4.21) 

The first sum contributes only if Aj, Aj belong to the same cluster, while the 
second contributes only if they belong to separate clusters. Because of the 
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geometrical mapping, ijj(Xi) = in the first sum. Also, in the sum over 

{A m : cl}, every cluster other than the one containing Aj and Aj contributes 
a factor s, the number of possible spin assignments. The cluster containing 
Xi and Xj , on the other hand, contributes a factor (s — l) 2 if Aj = Xj = 1 
and (— l) 2 for the (s — 1) other possible choices for Aj = Aj. Therefore, 

E i>(X l MX j M l ,j) = s N c- 1 [(s-l) 2 + (s-l)(-l) 2 ]n( i ,j) (4.22) 

{A m : cl.} 

where N c is the total number of cluster in the configuration. Hence, 

P 



Z p N(s - 1)N\ 



Jdl---dN E P(conn.) E HXi)ip(Xj)n(i,j) 



/ conn. \ {\ m -ci.} 
\ states / 



fdl-.-dN £ P^onnO^^-l + lM^j) 
^■^p^ J r conn. \ 

\ states; 



s JVc fi(A i ,A i )) (4.23) 



The second sum in the r.h.s of Eq. ( |4.21| ) contributes only if Aj, Aj belong 
to different clusters. The cluster containing Aj contributes a factor [s — 1) if 
Aj = 1, and a factor (—1) in all the other [s — 1) cases. The same holds for 
the cluster containing Aj. The N c — 2 remaining clusters contribute each a 
factor s. Hence, 

E i>(xMXj)[i-n(i,j)} 

{\ m :cl.} 

= s N °-* [{a - 1) + (s - 1)(-1)] 2 [1 - = (4.24) 

We can now substitute the results Eqs. ((123]) and ( pj~24l) into Eq. ( ]4TT0| ), then 
take again the limit s — > 1. Repeating the steps leading to Eq. (f4.18|) , we 
end up this time with 

l-S^l + ^y dxdy(^2^(hj)m-mr 3 -y)) (4.25) 



We define now the function 



g\x,y) = 1 (E^uXfi-^-jD) (4.26) 
p(x)p(y) \^ / p 
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to be the pair connectedness function [p(x) is the density at position x\. From 
the definition of the meaning of g^ is 

p(x)p(y) g\x, y) dxdy = Probability of finding two particles in regions 

dx and dy around the positions x and y, such 
that they both belong to the same cluster. 

(4.27) 

which shows the pair-connectedness to be a generalization of the correspond- 
ing function in lattice percolation |I6 |. Hence, 

S = 1 + ^J dxdy p{x)p{y)g\x,y) (4.28) 

Usually, the system is translationally invariant, so that g'(x,y) = g'(x — y) 
and p(x) = p(y) = p. Then 

S=l + ^p 2 J dfg\r) = l+p J dfg ] {r) (4.29) 

Comparing Eq. ( f4.28|) with Eq. ( |4.18|) , we finally obtain the important rela- 
tionship 

g j (x,y) = lim 1 mT^ Yl S i>{a)i>{l)p {2 \x, a; y, 7) (4.30) 

Let us now introduce the spin pair-correlation function, defined as 

gP(x,a;y,j) = . u 1 . . p (2) (x, a; y, 7) (4.31) 
which tends to 1 when \x — y\ — > 00. We can now rewrite Eq. ( |4.30 ) as 

g\x,y) = lim — - ^ ^ip(a)ip(i) 9?{x, a; y, 7) (4.32) 

or 



g j (x,y) = lim^-^^ 2 (a)^ 2) 



x,a;y,a) 

' ^ * 7 (2), 



+ lim ^(«)^(7)^ 2) (^>«;f ; 7) (4-33) 

s — 1 ~T 
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We now have 

J2^ 2 (a)gi 2 \x,a;y,a) = (s - lfgf\S, 1; y, 1) + (s - 1){-I) 2 gf\x , a- y, a) 

(4.34) 



where a is any value of the spin different from 1. As a result 
f 



lim— $> 2 (a)<7< 2) . 



lim g( 2) (x,a;y,a) 



(4.35) 



where a ^ 1. 
Similarly, 

X! ^(«)^(7)^ 2) (^«;y>7) = (s - 1) E^Wfl 2 '^- 1 ;^) 

+ ( S -l)^(a) S f)(i>;f,l) 

- X X ^(«)^(7) ^ 2) (^, «; y, 7) 

= (s - l) 2 \gf\x, 1; y, a) + #(f, a; £ 1)" 
+ ( S -l)( S -2)#(f,a;y,r ? ) (4.36) 

where a 7^ rj are any values of the spin which are both different from 1. As 
a result, 

lim — !— X t%)t%) 9 { s 2 \x, a; y, 7) = - lim </? } (^ ^5 £ ^) ( 4 - 37 ) 

8-+1 S — 1 ~ 8-*l 

Substituting Eqs. ( |4.35| ) and ( |4.37| ) into Eq. ( |4.33| ), we have, finally, that 

g ] (x, y) = lim [^ 2) (x, a; y, a) - g {2) (x, a; y, 77)] (p < p c ) (4.38) 

where a, r) 7^ 1 and a 7^ rj. 

This last equation is easily understood in terms of the geometrical map- 
ping. If two particles x and y belong to the same cluster, they must have the 
same spin, say a spin a. Hence g\x,y) must be contained in g^\x, a; y, a). 
However, g^ 2 \x, a;y, a) also includes the case where the two particles have 
the same spin but belong to different clusters. This happens if the two clus- 
ters have been assigned, by chance, the same overall spin. However, since 
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such spin assignment is random, the probability of the two particles having 
the spins a and a (identical) is exactly the same as their having spins a and 
T], where now o ^ 77 (note that such an assignment automatically necessitates 
that the particles belong to different clusters). Hence (x, a;y, a) exceeds 
g^(x,y) by precisely gj®(x, a;y,rj). This is the meaning of Eq. (|4.38| ). It 
allows us to calculate the pair-connectedness by working out the spin pair- 
correlations and then taking the limit s —*■ 1. This completes the relation 
between the percolative quantities and those of the Potts fluid. 

5 Conclusion 

This paper has focused on the formal basis of continuum percolation theory. 
It has provided a non-perturbative definition of the fundamental quantities 
of the theory as well as showing formally how the binding criterion p(r) and 
the interaction v(r) enter them. We have seen that the quantities of interest 
in continuum percolation can be obtained from the s — > 1 limit of the Potts 
fluid. Specifically, the magnetization and the susceptibility become in this 
limit the percolation probability and the mean cluster size, respectively, while 
the pair-connectedness is, in this limit, the difference between two Potts pair- 
correlation functions. 

The advantage of this mapping is that the Potts fluid has a Hamilto- 
nian formulation. This is the key to applying the techniques of equilibrium 
statistical mechanics and phase transitions to the problem of continuum per- 
colation, a task which will be undertaken in future papers in this series. The 
line of attack is always to perform all calculations or theoretical derivations 
within the Potts fluid model, then to take the limit s — > 1 and thus obtain 
the corresponding values or expressions in the percolation system. 

The first such technique that one would think to apply to any problem 
of phase transition is mean field theory. The Hamiltonian formulation of the 
Potts fluid allows a mean field approximation to be defined, while it would 
be far from obvious how to do this directly for the continuum percolation 
system. Mean field theory turns out to be non-trivial to derive because of 
the continuum nature of the system and the presence of interactions. It is 
the subject of the next paper. 
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